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ABSTRACT 

We carried out a search for M-type Vega-like stars by correlating the IRAS Faint Source Catalog 
with Hipparcos selected M-type stars. Three stars with apparent IRAS 25 /im excess emission are 
shown instead to be non-IR-excess stars from ground-based 11.7 and 17.9 ^m photometry. Two stars 
previously suggested to have Vega-like mid-IR excess are also shown to be non-excess stars. These 
results imply that other suggested mid-IR excess stars in the literature may also be false excess stars. 
Detection threshold bias is apparently responsible for these bogus IR excesses. Sixty micron excess 
emission from a previously known M-type Vega-like star (GJ 803) is identified again. 

Subject headings: (stars:) planetary systems: proroplanetary disks — (stars:) circumstellar matter — 
stars: late-type — stars: individual (GJ 803) 
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1. INTRODUCTION 

During the past fifteen years, about two dozen papers 
have been published that describe searches for stars with 
excess infrared (IR) emission (for recent reviews on these 
"Vega-like" stars, see Lagrange et al. 2000 and Zuckerman 
2001). These searches employed different techniques for 
cross correlating IR and stellar sources with no consistent 
definition of what defines an IR excess (see Song 2000 for a 
summary). To date, about 400 Vega-like stars have been 
identified. The Vega phenomenon overlaps with some very 
important solar system formation epochs; gas giant planet 
formation at < 10 Myr and terrestrial planet formation 
at ~ 100 Myr. Thus, knowing stellar ages of Vega-like 
stars is essential to studying extra solar planetary system 
formation in detail. Consequently, there have been many 
efforts to estimate the ages of Vega-like stars (see Song 
2000; Spangler et al. 2001, and references therein). Stel- 
lar ages can be fairly accurately and relatively easily deter- 
mined for late-type stars (e.g., Song et al. 2000), however, 
ages of early-type stars are less reliable (Song et al. 2001). 
Most of the currently known Vega-like stars are early-type 
because more luminous stars produce larger IR fluxes for 
stars with the same quantity of dust and the surveys con- 
duced by IRAS were flux limited. Thus, it is desirable to 
increase the known number of late-type Vega-like stars 
for more precise age estimates of such stars. Additional 
identification of late-type Vega-like stars is very useful 
to statistically strengthen studies such as planetary for- 
mation in different environments and dust lifetime as a 
function of stellar mass and luminosity (see Song's (2000) 
suggestion of dichotomy of Vega-like stars for example). 
Also, late-type Vega-like stars are excellent laboratories 
for studying the early evolution of our solar system. 

Despite the great number of M-type stars compared to 
earlier types, only two with IR excesses, GJ 803 (Tsikoudi 
1988) and Hen 3-600 (de La Reza et al. 1989; Jaywardhana 



et al. 1999) have been identified to date. Almost all pre- 
vious studies (exceptions include Aumann & Probst 1991 
and Odenwald 1986) searched only for far-infrared (60 /im) 
excess and were limited by the IRAS ~90% completeness 
sensitivity of 280 mJy in that band. IRAS sensitivity at 
25 /im was 210 mJy, and an MO star located 10 pc away 
with 200 K dust grains absorbing and emitting like black- 
bodies such that L/fl/Lstar=l-5 xlO^^, would have been 
detected by IRAS at 25 /im (Pstar+disk — 250 mJy) but 
not at 60 /zm {Ystar+disk - 100 mJy). 

A more complete survey for dust must be able to de- 
tect the stellar photosphere at high precision in order to 
evaluate excess emission above that level. One can do a 
fairly thorough disk excess assessment for stars whose pho- 
tospheres were detected by IRAS. However, around stars 
whose photospheres were too faint to be detected, only un- 
usually large disk excesses could be detected. At 60 /im, 
IRAS could detect the photosphere of an AO star out to 
~20 pc, but an MO only out to ^2 pc. The situation is 
somewhat better at 25 /im, where IRAS could detect the 
photosphere of an AO star out to ~50 pc and an MO pho- 
tosphere out to ~5 pc. By searching the catalog only at 
60 /im, previous surveys have not probed large regions of 
phase space where excess may exist around late-type stars. 

In this study, we concentrated mainly on M-type stars 
and attempted to perform the most thorough search for 
M-type Vega-like stars to date, especially at 25 /im, based 
on the IRAS FSC (Moshir et al. 1992) and Hipparcos cat- 
alog (Perryman et al. 1997). As a check on recently re- 
ported infrared excess stars, however, we also report on 
the F-type star HD 2381. 

2. SEARCH 

Based on the Hipparcos catalog, we selected ^ 530 
nearby (< 25 pc) stars with [B — V) + o-{b-v) > 1-40, 
where a^B_v) is uncertainty of {B — V). The Hippar- 
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cos catalog contains almost all early M-type (MO-2) stars 
within 10 pc. Many IR sources in the FSC with opti- 
cal stellar identifications are giant stars (Zuckerman et al. 
1995; Odenwald 1986). Therefore, one needs luminosity 
class information to identify Vega like stars. Following 
Silverstone (2000), we used a constraint on the absolute 
visual magnitude (My > 7.5 x {B — V) — 5.0) to ensure 
that a candidate is not a giant star whose IR-excess mech- 
anism may be different from that of a dwarf. Six stars from 
our initial sample do not meet the absolute visual magni- 
tude cut and they are HIP 21421, 50798, 66212, 66906, 
75187, and 82099. HIP 21421 and HIP 66212 are K-type 
giants and the other four stars appear to be main sequence 
stars with large uncertainties in {B — V). Among the four 
rejected main sequence stars, HIP 75187 is the only one 
detected by IRAS (only at 12 /um) and the measurement 
agrees with the flux density expected from its photosphere 
alone. Then, our sample stars were cross correlated with 
FSC sources with a maximum allowed offset of 30" be- 
tween Hipparcos and IRAS source positions (both at epoch 
1983.5 and equinox 1950). Only 152 stars from our initial 
sample have IR counterparts; among them, 96 stars were 
detected only at 12 /im and 55 stars were detected at both 
12 and 25 /xm. GJ 803 (AU Mic) was detected at 12 and 
60 /xm and was the only dwarf M-type star that has been 
detected at the IRAS 60 /um band. No objects were de- 
tected at 100 /xm. A previously known M-type IR excess 
star. Hen 3 600 (TWA 3), is not identified in our search 
because it is not bright enough to be included in the Hip- 
parcos catalog. 

Positions given in the published IRAS catalog are 
weighted means of 12, 25, 60, and 100 ^m source po- 
sitions based on their signal-to-noise ratios. Sometimes, 
25 and 60 fim sources are background objects far away 
from the stellar 12 /xm sources. Thus, in IR-excess sur- 
veys, it is mandatory to check each band's source position 
and to confirm that source positions in each band are co- 
incident. We checked the offsets among 12 and 25 fim 
positions (60 /xm position also for GJ 803) by using the 
"LONG FSC" from The Infrared Processing and Analysis 
Center (IPAC) at Caltech and foimd that only one object 
(GJ 433) shows a substantial (36") offset between the 12 
and 25 fim IRAS sources. For comparison, a median posi- 
tional uncertainty of IRAS FSC sources in the cross scan 
direction is 20". Offsets between the 12 and 25 /tm source 
positions for the other stars are negligible with respect to 
the IRAS positional uncertainties. 

To identify IR excesses, we performed spectral energy 
distribution (SED) fitting by using all known photometric 
data from the literature (queried through SIMBAD) in- 
cluding onhne 2MASS data. For GJ 413.1 and GJ 433, 
JHK magnitudes (see Table 1) were measured on 24 
November, 2000 (UT) with the NASA IRTF telescope at 
Mauna Kea Observatory. Since the accuracy of the pho- 
tospheric flux estimation at 12 /im depends strongly on 
the availability of near IR photometric data (i.e., JHK 
magnitudes), we have not carried out a SED fit to stars 
with only 12 /tm detections because most of these stars 
lack near IR photometry. Stellar SEDs are different from 
that of a blackbody. Opacity sources absorb light at wave- 
lengths with high opacities and re-radiate it at wavelengths 
with relatively low opacities. This results in a spectral en- 



ergy distribution very different from that of a black-body. 
Therefore we used PHEONIX NextGen synthetic stellar 
spectra (Hauschildt et al. 1999) instead of a blackbody 
SED. Among three SED fitting parameters (parallax, stel- 
lar radius, and effective temperature), Hipparcos parallax 
has been treated as constant. Stellar radius and effective 
temperature were estimated from the Hipparcos [B — V) 
value by using the spectral type versus colors/Te///radius 
relation of de Jager & Nieuwenhuijzen (1987). 

To quantify the strength of IR excesses, we defined r, 
the specific IR excess, ("specific excess" hereafter) as 

pIRAS _ pest 

25 /[im 25 /[im \ 

= pest V-"-; 

-'^^25 Aim 

where F^i^j^ and F^§*^^ are IRAS FSC 25 /xm flux and es- 
timated photospheric contribution at 25 /<m, respectively. 
For GJ 803, we used Fff(f^^,f and F|f^^ to calculate its 
60 jjm specific excess {reofim)- 

As shown in Figure 1, we found three stars (GJ 154, 
413.1, and 433) with r > 2.0 based on 25 /xm fluxes and 
a different star (GJ 803) with r = 7.60 based on 60 /tm 
flux. Contrary to an expected median specific excess value 
of zero for non-IR-excess stars, Figure 1 shows a median 
value of ~ 0.1 which may be due to the 25 /tm flux overesti- 
mation as explained in the IRAS Explanatory Supplement 
Version 2, III 131. AU IRAS flux density values in Table 1 
and Figures 2 3 are color corrected using Table VI. C. 6 of 
the IRAS Explanatory Supplement Version 2. For 12 and 
25 /tm fluxes, stellar effective temperatures were used to es- 
timate color correction factors. However, for 60 /im fluxes, 
if any IR excess exist (e.g., GJ 803), then dust tempera- 
tures were used instead of stellar effective temperatures. 

3. GROUND-BASED MID-IR PHOTOMETRY 

Mid-infrared imaging was performed with the facility 
instrument, the Long Wavelength Spectrograph (LWS) 
(Jones & Puetter 1993), on the 10 m Keck I telescope on 
UT 11 December 2000 and 4 5 February 2001. During aU 
three nights, the weather was photometric with low water 
vapor optical depth. LWS uses a 128 x 128 pixel Boeing 
Si:As detector, and has a plate scale of 0.08 arcsec/pixel, 
resulting in a focal-plane fleld of view of 10. 24 x 10. 24. 
Each object was measured in filters centered at 11.7 /tm 
(FWHM=1.0 /im) and 17.9 /xm (FWHM=2.0 /xm). Images 
were obtained at four positions by chopping the secondary 
at 2.5-5 Hz with a throw of 10" and nodding the telescope 
10" after ~ 20 s. In basic data reduction, the images 
were double differenced to remove the sky and telescope 
background, and bad pixels were corrected by interpola- 
tion. Throughout the nights, including just before and 
after each of the M-star measurements, bright infrared 
standard stars were observed for photometric calibration. 
Standard star measurements over the whole of each night 
were averaged and the standard deviation in their photom- 
etry was used as an estimate of calibration uncertainty. 
On 11 December, the uncertainty in the calibration was 
5% and 6% at 11.7 and 17.9 /tm, respectively. On 4 and 5 
February, the uncertainties were 15% at both wavelengths. 

Photometry was performed in a 16 pixel (l."3) diame- 
ter synthetic aperture on each image and the results are 
reported in Table 1. For an M-star of luminosity 0.1 Lq, 
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blackbody-like grains at a thermal equilibrium tempera- 
ture of 200 K will sit 0.6 AU from the star. Therefore, at 
a distance of 10 pc, any 12 or 18 /J.m excess should appear 
< 0. 1 in size, or spatially unresolved. It is clear from Fig- 
ure 2 that the apparent 25 /J,m IRAS excesses of GJ 154, 
413.1, and 433 are not real. We interpret this discord as 
follows: 

For the faint stars under consideration whose real fluxes 
are near detection threshold, a downward noise fluctua- 
tion could place the 25 ^m fluxes below the IRAS detec- 
tion threshold; thus none displays a significant 25 fxm flux 
deficit (negative r). Occasional large upward noise fluc- 
tuations could boost 25 fiia fluxes so that they would be 
classified as IR excess stars (positive r, "Detection thresh- 
old bias" or "Malmquist bias"). The final configuration 
thus resembles our Figure 1, with some excess stars but 
with no significant deficit star. In fact, the IRAS 25 fim 
S /N ratios of all of our three false IR excess stars are ~ 4 
which is the IRAS threshold value. 

4. STATISTICAL SIGNIFICANCE OF IR EXCESS 

Recently, Fajardo-Acosta et al. (1999, 2000) suggested 
that certain stars possess excess emission as measured by 
IRAS or ISO. We checked IR excesses at GJ 816 and 
HD 2381 with 11.7 and 17.9 /zm (18.7 /zm for GJ 816) Keck 
photometry. Apparent excesses for both stars turned out 
to be false positives (Figure 3). GJ 816 is not an IRAS 
FSC source and Fajardo Acosta et al. (1999) used Infrared 
Space Observatory (ISO) data. An incorrect ISO flux cal- 
ibration (for GJ 816) and Malmquist bias (for HD 2381) 
similar to our three false IR excess stars may be respon- 
sible for these apparent excesses. 

An occasional large upward noise fluctuation (e.g., 2a ~ 
2% probability) does not significantly infiuence stars with 
high signal-to-noise ratios; however, it can significantly af- 
fect stars with low signal-to-noise data. For our initial 152 
IRAS sources, we expect ~ 3 to have flux overestimates 
> 2a, in apparent agreement with what we have found. 
Based on this fact, some suggested Vega-like stars — gen- 
erally identified through huge surveys often encompassing 
thousands of input stars — could also be non-IR-excess 
stars. Thus, we suggest the following criteria for bona-fide 
Vega-like stars; (1) high S/N not subject to a Malmquist 
bias, (2) low S/N detections at 2 or more wavelengths. 



or (3) ground/spacc-bascd confirmation (e.g., Silverstone 
2000 and this study) with higher sensitivity and better 
spatial resolution than IRAS. 

5. SUMMARY AND DISCUSSION 

We have performed a search for IR excess emission 
among M-type stars by correlating the IRAS Faint Source 
Catalog with Hipparcos selected late type stars. Besides 
the previously known Vega-like star (GJ 803), three tenta- 
tive excess stars were identified, but these excesses turned 
out to be false based on our ground-based mid-IR pho- 
tomtery. Detection threshold bias (Malmquist bias) is 
thought to be responsible for these bogus IRAS IR ex- 
cesses. Two other stars (GJ 816 and HD 2381), suggested 
to be Vega like in the literature, are also shown to be non- 
IR-excess stars. In future studies, one should be aware 
that some Vega like; stars reported in the literature with 
low S/N ratios may be non-IR-excess stars as well. This 
is likely to be the case for most stars listed by Fajardo- 
Acosta et al. (2000). 

GJ 803 and Hen 3-600 show strong 60 /im excesses; 
they are the only unambiguously identified M-type dwarf 
stars with IR excesses. This could be due to the extreme 
youth of GJ 803 (12 Myr, Zuckerman et al. 2001) and 
Hen 3-600 (8 - 10 Myr, Webb et al. 1999). Song et al. 
(2001) have found two very young (~ 12 Myrs) late-type 
stars (HIP 23309, MO and HIP 29964, K6) co-moving with 
P Pictoris. Even if one assumed that these two Hipparcos 
stars have the same fractional IR luminosity as f3 Pictoris 
{Lju/Lstar ~ 10~^), their corresponding 60 /um fluxes 
(< 80 mJy and < 40 mJy, respectively) are below the 
IRAS detection threshold. This is true for late-type stars 
with (3 Pic-like excess in nearby young stellar groups, i.e, 
TWA. These stars would be excellent targets for future IR 
excess surveys by SOFIA or SIRTF. 

We are grateful to Mr. Michael Schwartz for assistance 
obtaining JHK data for GJ 413.1 and GJ 433 and Dr. M. 
Jura for helpful discussion and assistance with a Keck ob- 
servation. We also thank an anonymous referee for sugges- 
tions that clarifled the paper. This research was supported 
by the UCLA Center for Astrobiology and by NASA. We 
have used the SIMB AD/ Vizier database. 
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Table 1 
m-type ir-excess candidates 



GJ 


Sp. 
Type 


dist 
(pc) 


near 


IR data ( 


mag)' 


IRAS flux (mjy) 


r 

value 


Keck flux (mJy) 


Prediction 


(mJy)" 


excess? 


J 


H 


K 


12 jum 


25 /im 


60 


11.7 /im 


17.9 Mm 


11.7 Aim 


17.9 Mm 


154 


MO 


14.6 


6.67(3) 


6.03(5) 


5.85(5) 


150±24 


114±53 


< 199 


2.18 


159±13 


74±11 


157 


68 


NO 


413.1 


M2 


10.7 


7.23(2) 


6.55(2) 


6.23(2) 


141±21 


80±20 


< 130 


3.38 


177±27 


43±20 


147 


64 


NO 


433 


Ml. 5 


9.0 


6.46(2) 


5.95(2) 


5.67(2) 


205±25 


108±27 


< 101 


2.28 


213±11 


101±20 


214 


93 


NO 


803 


MO 


9.9 








537±32 


< 215 


273±46 


7.60+ 






574 


257 


YES 


816 


M3 


13.8 


7.55(1) 


6.96(2) 


6.69(2) 










93±13 


55^ 


101 


45 


NO 


HD 2381 


F2V 


74.2 


6.99(1) 


6.86(4) 


6.74(1) 


127±33 


< 73 


< 170 


2.30* 


55±5 


28±14 


55 


24 


NO 



* near IR data for GJ 154 from Alonso et al. (1994), for GJ 413.1 and GJ 433 from our IRTF measurements, and for GJ 816 and HD 2381 from 2MASS database. 
** expected photospheric flux. 

^ 60 /Ltm spec, excess. 25 /Ltm value is upper limit. 
12 fj,Tn spec, excess. 25 ^^m value is upper limit, 
this is 18.7 fim flux upper limit not 17.9 fim. 
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Fig. 1. Histogram of specific IR excess (r25/am.) for the 55 stars discussed in the text. The r-value of GJ 803 (reo/am = 7.6) is outside of 

the displayed range. 
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Fig. 2. — Spectral energy distribution fits of M-type stars with tentative IR excesses identified from this study by using the IRAS FSC. 
Solid circles are JHK and IRAS data, and diamonds indicate our ground— based 11.7 and 17.9 ^m fluxes. Thin solid lines are synthetic stellar 
spectra fit to visual and near IR (A < 2 fim) photometry {[M/H] = 0.0 and logg = 5.0) and a dotted line (only for GJ 803) indicates a dust 
component with T = SO K and Lir/L* = 6.7 X 10^*. Wavelength and flux density scales are logarithmic. Horizontal bars across JHK and 
IRAS data points indicate passband widths. 
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Fig. 3. — Spectral energy distribution fits of Vega-fike stars from the Hterature. Symbols have the same meaning as in Figure 2. For 
GJ 816, open squares show ISO fluxes (there is no IRAS data) and open diamonds are LWS measurements with the Keck I telescope. 



